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Abstract: The problem of decoupling a linear system by dy- 
namic compensation into multi-input multi-output subsystems is 
studied by applying proper and stable fractional representations of 
transfer matrices. A necessary and sufficient condition is given for 
a decoupling and a stabilizing controller to exist. The set of all 
controllers that decouple and simultaneously stabilize the system is 
determined in parametric form. Optimal decoupling controllers are 
then obtained by an appropriate selection of the parameter. 
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1. Introduction 


Decoupling is a way to decompose a complex system into non-interacting subsys- 
tems. In fact, certain applications necessitate controlling independently different 
parts of the system. Even if this is not required, the absence of interaction can 
significantly simplify the synthesis of the desired control laws. 

The decoupling problem has received much attention in the literature. For 
linear systems, different approaches have been used and control laws of various 
structure and complexity applied. 

The basic form of decoupling into single-input single-output subsystems is 
often referred to as the diagonal decoupling. This problem was posed by Vozne- 
senskij (1936) and studied by Kavanagh (1957), Strejc (1960), Mejerov (1965), 


*Submitted: May 2012; Accepted: October 2012 

+The original version of this paper was presented at the 18th International Conference on 
Process Control that was held in Tatranska Lomnica, Slovakia during June 2011. 

*Supported by the Technology Agency of the Czech Republic, Competence Centres Pro- 
gramme, project TE01020197 Centre for Applied Cybernetics 3 


140 V. KUCERA 


and Wolovich (1974). The studies were related to the inversion problem of ra- 
tional matrices. Attention was paid to the existence of proper rational transfer 
matrices. The issue of stability, however, was not properly addressed. 

A deeper insight was provided by the state-space approach. The pioneering 
work is due to Morgan (1964), who posed the problem of decoupling by static 
state feedback. Falb and Wolovich (1967) established a solvability condition 
while Gilbert (1969) related this condition to state feedback invariants of the 
system. Descusse and Dion (1982) then interpreted this condition in terms of 
system’s structure at infinity. 

The use of restricted static state feedback, namely the static output feedback, 
in decoupling was studied by Howze and Pearson (1970), Howze (1973), Denham 
(1973), Hazlerigg and Sinha (1978), Filev (1982b), Descusse and Malabre (1982), 
and Descusse, Lafay and Kuéera (1984). This is a very restricted problem, whose 
solution is hard to obtain, but it is very useful in applications. 

A more general form of decoupling into multi-input multi-output subsystems 
is referred to as the block decoupling. This problem was introduced by Wonham 
and Morse (1970) and Basile and Marro (1970). Using a geometric approach, 
they determined the solvability of the problem by static state feedback in several 
special cases. An alternative algebraic approach based on the structure algo- 
rithm was presented by Silverman and Payne (1971). Relationships between the 
two approaches were studied by Filev (1982a). 

The decoupling by dynamic state feedback was studied via the geometric 
approach by Morse and Wonham (1970), who obtained a deep insight into the 
internal structure of the decoupled system. By this time, the problem of decou- 
pling by dynamic state feedback was solved, including stability or pole distribu- 
tions that may be achieved while preserving a decoupled structure. The status 
of non-interacting control was reviewed by Morse and Wonham (1971). 

A comeback of the transfer function methods in the study of block decou- 
pling is witnessed through the works of Koussiouris (1979), Hautus and Hey- 
mann (1983), and Kuéera (1983). A dynamic state feedback was shown to 
be equivalent to combined dynamic output feedback and feedforward reference 
compensation, often referred to as a two-degree-of-freedom controller. To ad- 
dress stability issues, the Youla-Kucera parameterization of all stabilizing con- 
trollers was invoked. The basic results are reported by Kuéera (1983), Hautus 
and Heymann (1983), and Gémez and Goodwin (2000). The class of all de- 
coupled transfer matrices that can be achieved by a stabilizing controller was 
parameterized by Desoer and Gündeş (1986) and Lee and Bongiorno (1993). 
This result has made it possible to derive the H2-optimal decoupling controller, 
which minimizes the performance deterioration due to decoupling. 

The two-degree-of-freedom controller structure is ideally suited to decou- 
pling since only one of the degrees of freedom is affected by the decoupling 
requirement. This is not true for a pure feedback, or a one-degree-of-freedom 
controller. This case is considerably more difficult to solve, as shown by Hammer 
and Khargonekar (1984), Lin (1997), Youla and Bongiorno (2000), Bongiorno 
and Youla (2001), and Park (2008a). 
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Finally, the decoupling in the generalized plant model, which covers a broad 
range of control problems in a unified setting, was considered by Park (2008b). 
Such a plant model can accommodate non-square plant and non-unity feedback 
cases with one-degree-of-freedom or two-degree-of-freedom controller configura- 
tion. The benefits of such a general problem formulation consist in a unified 
treatment rather than in simplicity of the solution. Indeed, matrix operations 
need to be converted to vector operations with vectors of a much larger dimen- 
sion, which result from the Kronecker and Khatri-Rao products of matrices. 

This paper adopts the most general setting that is meaningful for decoupling: 
a system in which the measurement output may be different from the output 
to be decoupled and a dynamic controller that features both feedback and feed- 
forward parts. The class of all such controllers that decouple and stabilize the 
system is determined in parametric form and the parameter is used to obtain 
the Hz-optimal controller. The solution is simple and direct. The controller 
configuration implies that decoupling and stability are two independent issues. 


Figure 1. Control system 


2. Problem formulation 


Consider a linear, time-invariant, differential system governed by the input- 
output relation 


y= Syu, (1) 


where u is the q-vector input, y is the p-vector output and Sy is the transfer 
matrix of the system. We assume that Sy is a proper rational matrix over R(s), 
the field of rational functions. 

Let pı, ..., pk be a given set of positive integers that satisfy 


k 


a Pi =P: 


System (1) is said to be decoupled, or more specifically (p1, ..., p~)-decoupled, if 
there exist positive integers q1, ..., qk Satisfying 
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such that Sy has the block diagonal form 


Sy 


Sk 


where S; is pi X qi- 

This is not a generic property of the system, but it can be achieved by a 
suitable compensation. To this effect, let z denote the m-vector output ofthe 
system that is available for measurement and let it be related with the input by 
the equation 


z= Szu, (2) 


where S, is a proper rational matrix over R(s). 
The most suitable linear, time-invariant, differential controller can then be 
described by the equation 


u = K,v+ Kzz, (3) 


where v is an external reference input of appropriate dimension, say r. As it is 
seen in Fig. 1, the transfer matrices K, and K, represent the feedforward and 
the feedback parts of the controller, respectively. We assume that both K, and 
K, are proper rational matrices over R(s). 

The decoupling problem is then to find matrices K, and K, such that the 
transfer matrix 


T = S&I — K,Sz)"'Ky (4) 


from v to y be suitably block diagonal. 


Figure 2. Control system with the complete set of independent inputs and 
outputs 
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Obviously, unless additional provisions are made, the decoupling problem 
is trivial as it could be solved by K, = 0. Thus, it is necessary to impose 
certain admissibility condition on the decoupling controller to make the problem 
meaningful, for example 


rank T= rank Sy (5) 


over R(s). This condition is equivalent to the preservation of the class of con- 
trolled output trajectories. We thus require that no essential loss of control 
occurs through the decoupling process. 

Another requirement, frequently imposed on the decoupled system in prac- 
tice, is that of stability. This requirement means that the states of the system 
go to zero from any initial value. 


3. Preliminaries 


A stable system gives rise to a proper and stable transfer function. In order 
to study stability of the decoupled system it is convenient to express the trans- 
fer matrices of the given system and those of the controller in the following 
factorized form 


Is |=[e]a 


are proper and stable rational matrices that are right coprime and 
P,| -Q R] 


are proper and stable rational matrices that are left coprime. 

These proper and stable fractional representations exist and are unique up 
to right and left multiplication, respectively, by a unimodular matrix. Recall 
that a proper and stable rational matrix is said to be unimodular if its inverse 
exists and is proper and stable. 

The system equations (1) and (2) and the controller equation (3) then take 
the form 


J-J o 


u=P | -Q EE (7) 


v 
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The overall system transfer function reads 
T=C(PA+QB)"'R. (8) 


The fundamental assumption we make here is that the part of the given system 
that is not controllable from u is stable and the part of the given system that 
is not jointly observable from y, z is stable. Similarly, we assume that the 
controller is realized in such a manner that its part that is not jointly controllable 
from v, z is stable and its part that is not observable from u is stable. 

The issue of stability of the overall system is then solved as follows. 

Lemma 1. The overall system described by (6) and (7) is stable if and only 
if the matriz PA + QB is unimodular. 

Proof. In the overall system, inject inputs x and w as shown in Fig. 2. Then 
the overall system is stable if and only if the nine transfer matrices between the 
inputs v, w, x and the outputs u, y, z, given by 


u A x 
z |= |B | (PA+QB)"[ P -Q R] | w 
y C v 


are all well defined and proper and stable rational. This statement follows from 
the assumption of stability of the uncontrollable and unobservable parts of the 
system. 

Now, in view of the coprimeness assumptions on A, B, Cand P, Q, R these 
transfer matrices are well defined and stable if and only if PA + QB is a unimo- 
dular matrix. O 


4. Problem solvability 


A necessary and sufficient condition will now be established for a system to be 
decoupled and stable. 
Based on the partition (pi, ..., pe), write 


Cy 
C=]: |, (9) 
Cr 


where C; is a p; x q submatrix. 

Theorem 1. Given system (1), (2) in fractional form (6) and partition (9), 
there exists an admissible controller (3) such that the overall system is 

(i) stable if and only if A and B are right coprime, 

(ii) decoupled if and only if 


k 
Se rank C; = rank C. (10) 
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Proof. (i) Let the overall system be stable. By Lemma 1, the matrix PA + QB 
is unimodular whence A and B must be right coprime. 

Conversely, let the matrices A and B of (6) be right coprime. Then there 
exist proper and stable rational matrices P and Q such that 


PA+QB=I1 (11) 


with P invertible and the inverse of P proper. 

Then, controller (3) in fractional form (7) that is defined by the matrices 
P and Q from (11) and by an arbitrary proper and stable rational matrix R 
satisfying rank CR = rank C is admissible since, by (8), 


rank T = rank CR = rank C = rank Sy. 


The resulting system (1), (2) and (3) is stable in view of Lemma 1 and identity 
(11). 
(ii) Let (7) be an admissible decoupling controller for system (6). Denote 


K := (PA + QB)" R. 
The block diagonal property of the matrix T then implies 
k 
rank CK = ae rank CiK 


and the admissibility of the controller gives 


rank CK = rank C;, 7=1, k. 


sees 


Therefore (10) holds. 
The sufficiency will be proved by constructing a suitable R. Denote 


ri := rank Cj, i=1,..., k. 


Then there exists a p; xX p; unimodular proper and stable rational matrix U; 
such that 


aua |, (12) 


where the rows of C} are linearly independent over R(s) and where the zero 


matrix has p; — ri rows and may be empty. If (10) holds, then 


C:= |: 
Ct 
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has linearly independent rows over R(s). Hence, there exists a q x q unimodular 
proper and stable rational matrix U’such that 


Dı 0 
C'U' = oe, : ; (13) 
Dy, 0 


where D; is an r; x r; diagonal proper and stable rational matrix and where the 
zero matrices have q — r columns with r defined by 


k 
ris > Ti. 
i=1 * 


Define an admissible controller (7) by the matrices P and Q from (11) and by 
the matrix R formed by the first r columns of U’. The transfer matrix (8) 


v, Ia | 


T=CR= a 2 (14) 

Ux Dk 
[o] 
is block diagonal. The resulting system is therefore decoupled and the external 
reference input v has dimension r. O 
The interpretation of these solvability conditions is as follows. Condition 
(11) corresponds to the stability of the subsystem of the given system that is 
not observable at the measured output z. Condition (11) calls for the linear 
independence of any two outputs of the given system that belong to different 
blocks. The solvability of the decoupling problem thus strongly depends on the 


partition (p1, ..., Pk), that is to say, upon the allocation of the outputs into the 
blocks. 


5. Controller parameterization 


When a decoupling and stabilizing controller exists, we shall parameterize the 
class of all such controllers. 

The control system (6), (7) is stable if and only if PA + QB is a unimodular 
matrix by Lemma 1. Thus, stabilization involves only the feedback part K, of 
the controller, which surrounds the measurement subsystem S,. As a result, 
the parameterization of K, amounts to the well-known Youla-Kucera param- 
eterization of feedback stabilizing controllers. For details, see Kuéera (1975), 
Youla, Jabr and Bongiorno (1976), Kuéera (1979), Desoer et al. (1980), and 
Vidyasagar (1985). 

Let P, Q be any solution pair of equation (11). Then the solution class of 
(11) is given by 


P=P+WB, Q=Q-WA, (15) 


Optimal decoupling controllers revisited 147 


where A and B are left coprime, proper and stable rational matrices such that 
A'B=BA" (16) 


and W is an arbitrary proper and stable rational matrix parameter. 
The class of all stabilizing proper rational K, is then obtained in the form 


K,=—-P71Q=-(P+WB)"'(Q-WA), (17) 


where the parameter W is constrained so that the inverse of P+W B exists and 
is proper rational. 

Once the control system (6) and (7) is stabilized, it is decoupled if and only 
if T = CR by (8). Thus, decoupling involves only the feedforward part K, of 
the controller. 

Partition the q x q unimodular matrix U'defined in (13) as 


U=[ UL UL], 


where U; has r columns and Uj_, has q — r columns and may be empty. The 
class of all decoupling proper rational K, is then given by K, = P~'R with P 
determined in (14) and 


vi 
R=U, f (18) 


Vk 


where V; is an arbitrary r; xr; proper and stable rational matrix parameter. The 
matrices V1, ..., Vk, in turn, parameterize the class of achievable block-diagonal 
transfer matrices (8) as follows 


Di 
Uı 0 Vi 


T= on 2a SR . (19) 
i A 


The parameterization of decoupling stabilizing controllers reveals that decou- 
pling and stabilization are two independent issues. That is why the controller 
described by (3) is called the two-degree-of-freedom controller. However, this is 
no longer true for one-degree-of-freedom controllers, e.g., for the error-actuated 
controllers described by u = —P~'Q(v — w) in place of (7). 


6. Optimal controllers 


The decoupling constraint can deteriorate system’s performance. The bonus of 
having a parameterized solution set is that the lost performance can easily be 
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optimized. Optimal decoupling controllers can be obtained by an appropriate 
choice of the parameters V1, ..., Vk and W. 

Suppose that the control objective is for each block of outputs y; to track 
the corresponding block of reference inputs v;. Thus, we suppose that p; = ri 
for i = 1, ..., k, i.e., there are as many reference inputs as controlled outputs in 
each block. The tracking error for each block is given by 


€i = Vi — Yi = Hivi. 
In view of (19), H; has the generic form 


where F; := U; D; and V; are proper and stable rational matrices with F; fixed 
and V; an arbitrary parameter to be specified. 

The benefits of controller parameterization will now be demonstrated in 
the case of Hə control design, see Kučera (2011). It turns out that only the 
parameters Vi, ..., Vk are subject to selection, whereas W is free and can be 
independently selected to accommodate additional design specifications. 

Suppose that for each block, the reference-to-error transfer function H; is to 
have the least Hz norm defined by 


oo Ip 
| Hill, = (trace x | H; Go) Ha (jw) d 


where the asterisk denotes the conjugate transpose. Thus, H;(s) := HT (—s8) 
for any complex argument s. 
To achieve this task, determine the inner-outer factorization of F;, 


F; = Fir Fio, 


where Fj; is inner and Fio is outer. Since F; is square and nonsingular, Fj; 
satisfies F Fiz = Land Fio is free of zeros in Res >0. 
Since Fj; is inner, left multiplication by Fy preserves the H 2 norm, 


Hille = Fir Hilly = Fir — FoVilla - 


Observe that F; (00) = I. Separate the strictly proper part, F; 


* =] 
‘Tey’ of Fi as 
follows 


Fi =I+ Fitsp 
and note that, by definition, Fitsp has poles only in Res >0. Then 
2 x 2 
lE = || Ets + O- Fio vol 
* 2 2 
= || Farepll; E Il he FioVill3 
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because the cross terms contribute nothing to the norm. This is a complete 
square in which only the second term depends on V;. Therefore, a unique V; 
that attains the minimum of the norm for subsystem 7 is 


V; = Fo- (21) 


However, only a proper and stable V; is admissible. It follows that the H2 control 
problem for subsystem 7 has a solution if and only if Fio is unimodular. The 
minimum norm is then given by 


min || Hill, = ||Firspll2 - 


Other performance requirements can be addressed by the choice of W; for ex- 
ample, Kuéera (2012) studies asymptotic tracking of a reference. 


7. An example 
Consider a system defined by (1), (2) with the transfer matrices 


1 s+2 2s+1 3s 
= -—1 = 2 —1 
Sy a | s—1 E ’ Sz c | aur 5 f 
s+2 s+2 


Thus, the measurement output z is different from the output y to be decoupled 
in that it involves a non-unity feedback sensor. 

The task is to determine a two-degree-of-freedom controller (3) that (1, 1)- 
decouples and stabilizes the system. 

The first step is to obtain a proper and stable fractional representation (6) 
for the system. Standard calculations yield 


2s+1 3s 1 1 
-E B Lbs bal 

s+2 s+2 s+2 s+2 
Now apply Theorem 1. Since A is right coprime to B, a stabilizing controller 
exists. Since the rank of C equals the sum of the ranks of the rows of C, an 
admissible decoupling controller exists as well. 

All stabilizing and decoupling controllers will be parameterized using the 

fractional representation (7). To obtain the feedback part of the controller, we 
consider any particular solution of equation (11), for example 


S 1 0 = 0 0 
P= s , a=] le 
= al 1 0 


150 V. KUCERA 


A left coprime fractional representation that satisfies (15) is given by 


As | ia se B= 


0 


s—l 

s+2 2 

_ (s-1)(2s+1) 38 : 
s+2 s+2 


Thus, the solution class (15) of equation (11) is 


s—1 
1 0 s 2 
P= | _2stl _9 | +W | E-241) 35 | 
s+2 ~~ (s-£2)2 sF? 


a[i 8)" i] 


z s+2 


(22) 


To obtain the feedforward part of the controller, note that U; = U2 = 1 and the 
unimodular matrix defined in (13) equals 


,„ f2 -= 
palaa h 


Thus, (17) yields 


refa i) [0 &) (2) 


The matrices P, Q in (22) and R in (23) define the class of all controllers that 
solve the given problem. The parameters Vi, V2 are free proper and stable 
rational functions and W is permitted to range over proper and stable rational 
2 x 2 matrices so that the inverse of P exists and is proper. Obviously, this 
means that P(oo) is to be a nonsingular matrix. 

The decoupled transfer matrices that can be achieved in this example are 
given by (19) as 


1 0 v o0 
Pale callow] 
| 0 2 0 V 
The optimal controller that minimizes the Hə norm of the reference-to-error 


transfer matrix is determined from (21), channel by channel. Clearly, V; = 1. 
To optimize V2, the inner-outer factorization of 


and the strictly proper part of 
k s+1 


at s—1 
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equals 
2! 62 
2Isp — E 4 
Thus, from (21), 


s+2 


V = : 
i s+1 


It follows from (23) that the unique optimal R is 


2 —s+2 
[ha 


s+1 


and the overall system has the transfer function 


D 
s+1 


The parameter W does not affect the transfer function T. It is related to the 
internal structure of the decoupling controller. For example, W = 0 results in a 
fairly simple feedback controller. Another choice of W may provide a robustly 
stabilizing controller. 


8. Conclusion 


An optimal Hə decoupling control problem has been studied in the most general 
setting, for systems in which the measurement output may be different from the 
output to be decoupled and for dynamic controllers that feature both feedback 
and feedforward parts. The class of all such controllers that decouple and sta- 
bilize the system has been determined in parametric form and the parameter 
has been used to obtain the Hə-optimal controller. 

The main contribution of the present paper is in a streamlined and trans- 
parent exposition and a simple and direct solution of the optimal decoupling 
control problem, see Kučera (2011). This is primarily because of the following 
facts. The adopted controller configuration is ideally suited to decoupling since 
stability and noninteraction can be treated as two independent constraints. The 
problem is formulated and solved using an algebraic approach, namely the no- 
tion of proper and stable fractional representations for system transfer matrices. 
The parameterization of the decoupling controllers is achieved via the Youla- 
Kučera parameterization of all stabilizing controllers. Finally, the Hə norm 
involved in the optimization is minimized using the completion of the squares, 
which is a simple algebraic technique. 

A large body of literature exists on decoupling and related topics. In tech- 
nical details, the present paper draws inspiration from the work of Hautus and 
Heymann (1983) for the formulation of the problem, from Kučera (1983) for the 
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algebraic treatment of stability, from Desoer and Gündeş (1986) for the param- 
eterization of the decoupled system, and from Lee and Bongiorno (1993) for the 
optimal control. 

The matrix fraction approach presented in this paper can be also applied to a 
class of linear, time-invariant, singular systems that possess a transfer function. 
This is a subject of current research. 
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